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ABSTRACT 

We examine the properties of the outflowing matter from an advective accretion 
disc around a spinning black hole. During accretion, rotating matter experiences cen¬ 
trifugal pressure supported shock transition that effectively produces a virtual barrier 
around the black hole in the form of post-shock corona (hereafter, PSC). Due to shock 
compression, PSC becomes hot and dense that eventually deflects a part of the in¬ 
flowing matter as bipolar outflows because of the presence of extra thermal gradient 
force. In our approach, we study the outflow properties in terms of the inflow param¬ 
eters, namely specihc energy {£) and specific angular momentum (A) considering the 
realistic outflow geometry around the rotating black holes. We find that spin of the 
black hole (a^) plays an important role in deciding the outflow rate Rm (ratio of mass 
flux of outflow and inflow), in particular, is directly correlated with for the 
same set of inflow parameters. It is found that a large range of the inflow parameters 
allows global accretion-ejection solutions and the effective area of the parameter space 
{£, A) with and without outflow decreases with black hole spin (a^). We compute 
the maximum outflow rate as function of black hole spin (a^) and observe 

that R^°'^ weakly depends on ak that lies in the range ~ 10 % — 18% of the inflow 
rate for the adiabatic index ( 7 ) with 1.5 ^ 7 ^ 4/3. We present the observational 
implication of our approach while studying the steady/persistent Jet activities based 
on the accretion states of black holes. We discuss that our formalism seems to have 
the potential to explain the observed Jet kinetic power for several Galactic Black Hole 
sources (GBHs) and Active Galactic Nuclei (AGNs). 

Key words: accretion, accretion disc - black hole physics - shock waves - ISM: jets 
and outflows -X-rays: binaries. 


1 INTRODUCTION 


Powerful Jets and outflows are commonly observed in ac¬ 
creting black hole systems including activ e galactic nu¬ 
clei (AGNs) and X-ray binaries (XRBsi llMirabel et al\ 


I 992 I: iMirabel & Rodrieuedll994 iHiellming & Ruoenl 

1995; 

Ferrari 

1998: Mirabel & Rodriauez 19981: Junor et al\ 

199 c; 

Cheuna 

2 OO 2 I: Mirabelll2003l: Miller et all 20121). In snite of 


the availability of wealth of high resolution observations, 
the physical mechanism of Jet generation and its power¬ 
ing processes are still remain unclear. However, the obvious 
cause seems to be the extreme gravity that powers the out¬ 
flows where spin of the central objects may play an impor¬ 
tant role in producing the relativistic Jets. Earlier work of 
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iPenrosel (1 19691 ) demonstrated that infalling particles on to 
a rotating black hole have the potential to extract some of 
the rotational energy of the black holes. Close to the hori¬ 
zon, space-time geometry is dragged due to black hole ro¬ 
tation and the magnetic field lines are twisted resulting the 
transport of energy f r om th e black hole along the field lines. 
iBlandford &: ZnaiekI (Il97i1 ) showed that this purely electro¬ 
magnetic energy extraction mechanism via the threading of 
magnetic field lines has the potential to power the Jets. This 
appealing mechanism clearly indicates that black hole spin 
perhaps play an important role to produce Jets in black 
hole systems. In order to establish this feature observation- 
ally, attempts have been made to find the evidence of spin- 
po wered Jets fo r the stellar mass black holes . Acco rding 
to ISteiner et al. I ll2013l) and iMcClintock et all (|2014| ). sig¬ 
nificant positive correlation between the radio luminosity 
associated with the mechanical power of ballistic Jets and 
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the spin of the black hole is seen when the transient black 
hole syst ems are fed at Eddin gton- limited accretion rates . 
However, iRnssell et al\ (|201lj i and iFender fc Galld (|2014l ~l 
ruled out the direct evidence of any such correlation be¬ 
tween the transient Jet po wer and the black hol e spin. Using 
radio-optical correlation, Ivan Velzen &: Falckd (l2013l l also 
reported that the black hole spin possibly does not have any 
dominant role in poweri ng the Jets. In a numerical endeavor, 
iDe Villiers et oZl ll2005l ) claimed that the Jet efficiency in¬ 
creases by two orders of magnitude for rapidly rotating black 
hole co mpared to the non - spinn ing black hole. In the similar 
context. iFernandez et all (l2015ll showed that the more mass 
is evacuated from the disc as the black hole s pin is increased. 
In an early effort . [Ponea fc BiermannI lll996li also found that 
the Jet power is strongly dependent on the spin of the black 
holes without considering the Jet geometry. Overall, all the 
above findings are in contrast and therefore, inconclusive. 


Meanwhile, some authors pointed out that there exists a 
confirm connection between the accretion spectral states and 
the Jet states (i.e., thick an d thin) in black hole systems, par¬ 
ticula rly for microquasars llGallo et al.ll2003 : iRushton et al\ 
l2010l 'l. In the low-hard spectral states (LHS) of galactic black 
holes, quasi-steady and persistent Jets are observed whereas 
relatively stronger Jets are launched generally in the hard- 
intermediate state (HIMff). In addition, transient relativistic 
ejections are observed during the state transiti on from hard- 
intermedia te to soft-intermediate state (SIMS) llFender et al\ 
l2004l2009lh Interestingly, Jet activities are not seen in the 
high-soft states (HSS). It has also been observed that dur¬ 
ing ejections, Quasi-Periodic Oscillations (QPOs) are not 
observed and energy spectra get softens (IVad awale e t al\ 


_ spectra ge. _ - _- - 

I2OOII : iMiller-Jones et oZ.ll2012l : lRadhika fc Nandill2014| j. This 

is conjectured in a way that the matter is ejected from the in- 
ner part of the disc ( i.e., Compton coron a equivalently PSC) 

( Feroci et al\ 1 19991 : iNandi et al\ l2001al : Chakrabarti et m 
I2OO2I I. Towards this, very recently, Radhika fc Nandi ( 2014l l 
reported the inflow-outflow connection through the ob¬ 
servation in Radio and X-ray bands for the object XTE 
J1859-I-226 and claimed that the fast variation of the 
‘spectro-temporal’ properties of outbursting GBHs are con¬ 
nected with the disruption of the post-shock disc (i.e., 
PSC) in terms of ev acuation of matter in the form of Jets 
llRadhika et oZ.llioiSl ). They pointed out that the QPOs are 
not seen during ‘transient’ ejection of Jets observed in Ra¬ 
dio. This result indicates that P SC possibly be responsible 
for the origin of QPOs in GBHs llChakrabarti fc ManickamI 
I2OOOI : iNaiidi et aDl20oI3 f9). When Jets are emerged out 
at the cost of the evacuation/d i sruption of PSC, QPO s 
are missing (IVadawale et a/.ll200ll : iRadhika fc Nandlll20l4l . 
Extensive MHD simulations of accretion disk around ro¬ 
tating black holes also indicate that windy hot mate¬ 
rials (i.e., coron a) blows away fr o m the inner part of 


the disc as Jets ( Koide_eE^ I I2OO2I : iMcKinnev fc Gammid 
I2OO4I : [Pe Villiers et al\ 20051 ). Meanwhile, Das et al. ( 20141 ) 
showed that PSC modulates quasi-periodically when viscos¬ 
ity parameter is chosen to its critical value and such modu¬ 
lation successfully exhibits quasi-periodic variation of PSC 
as well as the outflow rates. These scenarios indicate that 
there seems to be direct correlation between PSC and out¬ 
flow. Based on these correlations, one can conjecture that 
the disk-jet symbiosis is strongly coupled and advective ac¬ 


creting disc perhaps be responsible for the launching of Jets 
and outflows. 


Advective accretion flow around the black holes must 
be transonic in order to satisfy the inner boundary con¬ 
ditions. Inflowing matter experiences a virtual barrier in 
the vicinity of the black hole due to centrifugal repulsion 
against gravity. According to the second law of thermo¬ 
dynamics, such a virtual barrier triggers the discontinu¬ 
ous transition of the flow variables in the form of shock 


possible (iFukud 1987: 

Chakrabartil 19891: Lu et al. 

I1999I: 

iBecker & Kazanad 2001 

: iFukumura & Tsurutal 12004 

). Due 


to compression, the shock induced accretion flow produces 
hot and dense post-shock corona (PSC) surrounding the 
black holes which essentially acts as the effective bound¬ 
ary layer of the black holes. During accretion, a part of 
the inflowing matter is deflected by PSC which is fur¬ 
ther driven out in the vertical direction by the excess 
thermal gradient force across the shock, producing bipo¬ 
lar outflows. Such an appealing mechanism to launch out- 


firmed through numerical simulations llMolteni et al\ 1994. 

199e 

: Machida et al. 

2OOOI: De Villiers et al\ 20051: Okudal 

2014 

:lDas et al.ll2014 

Okuda & Dadl2015l). In addition, nu- 


merous attempts have been made theoretica lly to calcu¬ 
late th e mass outflow rate around black holes. IChakrabartil 
l|l999l) calculated the mass outflow rate considering isother¬ 
mal fl ow. Using both Ke plerian and sub-Keplerian compo¬ 
nents, IPas et (l2001al ) self-consistently estimated mass 
loss from the disc and foun d that outflow rate strongly 
depends on both components. ISingh fc Chakrabartil (1201 ll ) 
implemented the energy dissipation across the shock while 
obtaining the outflow rate and found that possibility of 
mass loss anti-correla tes with the di ssipat ion rate. Fol¬ 
lowing the work of Molteni_et_^ ([19^), several au- 
thors JChntto2adh^a^_&_^a^_l2007UDas&_ChnttoMdh^ 


200^ Kumar fc ChattopadhvavT 20131 : Kumar et al\ 20131 : 


i 


Das et al.ll2014l) studied the properties of mass outflow rate 


in terms of inflow parameters considering dissipative accre¬ 
tion flow. However, all these works were carried out with 
limitation as spin of the black hole was not considered. 


Motivating with this, we model the inflow-outflow ac¬ 
tivity around a rotating black hole. Since the PSC is induced 
by the shock transition itself, we self-consistently calculate 
the properties of shock driven outflows in terms of the inflow 
parameters and investigate the effect of the black hole spin 
on it. For si mplicity, we adopt the pseud o-Kerr potential 
proposed bv IChakrabarti fc Mond^ (l2006l) to describe the 
space time geometry around the black hole. This potential 
accurately reproduces the particle trajectories around rotat¬ 
ing black holes having spin ak < 0.8. We calculate the global 
accretion solutions in presence of thermally driven outflows 
and obtain the parameter spaces spanned by energy and 
angular momentum of the inflowing matter in terms of the 
black hole spin. We observe that the shock induced global 
inflow-outflow solutions exist for wide range of inflow pa¬ 
rameters. Varying the inflow parameters, we estimate the 
maximum mass outflow rate as function of black hole 

spin (ofe) and find that R^^ weakly depends on a*, having 
highest value ~ 17% — 18% of the inflow rate for 7 = 4/3. 
Further, we calculate the unabsorbed X-ray flux of several 
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black hole objects using the RXTE^ satellite archival data 
and obtain their accretion rate s considering the accretion 
efficiency rj — 0.3 (ffh ornelllQT^ l. which perhaps reasonable 
for rotating black holes. With this, using our theoretical es¬ 
timate of maximum mass outflow rate, we then compute the 
maximum Jet kinetic power (1/™“). These results are fur¬ 
ther compared with those available from observations and 
close agreements are seen. Following this, our theoretical 
prediction essentially provides an estimate of for those 
sources for which observed Jet power is uncertain. 

In the next Section, we describe the basic assumptions 
and the governing equations for our model. In §3, we discuss 
the methodology to calculate the mass outflow rates self- 
consistently. In §4, we present the results which is followed 
by discussions. In §5, we apply our formalism to calculate the 
Jet kinetic power for several astrophysical objects (GBHs 
and AGNs). Finally in § 6 , we draw the concluding remarks. 


radial velocity and a is the adiabatic sound speed defined 
as a = Here, P is the isotropic pressure, p is the 

gas density, 7 is the adiabatic index, respectively. The ac¬ 
creting gas is described by the adiabatic equation of state as 
P = Kp'^, where K is the measure of specific entropy, which 
is constant except at the shock transition. The effe ctive po¬ 
tential 4? is given by dGhakrabarti fc MondallboOfil '). 


where. 


$ = - 


B + y/B^ - 4.AC 
2A 


A 


\2 
a A 


B^-l + 






-b 


2afeA 

r^x 


( 2 ) 


2 MODEL EQUATIONS AND ASSUMPTIONS 

We consider an axisymmetric disc-jet system around a ro¬ 
tating black hole in the steady state. Here, the accretion disc 
lies along the black hole equatorial plane while the Jet geom¬ 
etry is described about the black hole rotation axis. In this 
work, we mainly focus on the inner region of the accretion 
disc, where the viscous time scale is larger than the in-fall 
time scale and therefore, the angular momentum transport 
due to the differential motion is weakly significa nt leaving 
the flow to be inviscid there (IChakrabartil 1989l l. Further, 
in our model. Jets are originated from the inner part of the 
accretion disc (PSG) with the same angular momentum of 
the disc as we neglect the effect of resulting torque in the 
disc-jet system. Next, we present the governing equations 
that describe the fluid properties of the accretion disc and 
the Jets. All the equations are written in geometric unit sys¬ 
tem as G = Mbh = c = 1, where, G is the Gravitational 
constant, Mbh is the black hole mass and c is the speed of 
light, respectively. In this system, the unit of length, mass 
and time are expressed as GMbh /(?, Mbh and GMbh /(? ■ 


1 , 2 afe^ , 

0 — 1-^-^-—^—. 

r — xo X 2x^ 

Here, x and r represent the cylindrical and spherical radial 
distance considering the black hole is located at the origin 
of the coordinate system and A is the specific angular mo¬ 
mentum of the flow. Here, xo = 0.04 + 0.97afc -b O.OShal, 
uj = 2akl{x^ -b a\x + 2,a\) and = {x^ — 2x + a1)j{x'^ + 
Of. -b 2aflx), a is the redshift factor and represents 
the black hole rotation parameter defined as the specific 
spin angular momentum of t he black hole. According to 
IChakrabarti Sz Mondall (l2006l l. the above potential mimic 
the Kerr geometry quite satisfactorily for a wide range of 
fflfc < 0.8. 

(ii) the mass conservation equation: 

M = Anpvxh, (3) 

where, M denotes the mass accretion rate which is constant 
everywhere except the region of mass loss and h is the half- 
thickness of the dis k obtained from thin disk approximation 
llChakrabarti|[l989l l as. 


2.1 Governing equations for Accretion 

We consider a geometrically thin, axisymmetric, low angu¬ 
lar momentum, advective accretion flow around a rotating 
black hole. For simplicity, w e adopt the pseudo-Ker r effec - 
tive potential introduced bv iGhakrabarti k, Mond^ (l2006l f 
to represent the space-time geometry around the black hole 
instead of using full general relativistic prescription. This 
enables us to solve the problem following the Newtonian 
approach and keeping all the salient features of space-time 
geometry around it. The equation of motion describing the 
accreting matter is given by, 

(i) the energy conservation equation: 


S = 



7- 1 




( 1 ) 


where, S represent the specific energy of the flow, v is the 


^ http://www.heasarc.gsfc.nasa.gov 


where, = (cl4?/clr)z<<a;, 2 : is the vertical height in the 
cylindrical coordinate system and r = y/x'^ + z'^. 

Combining the expression of sound speed and the adi¬ 
abatic equation of the st ate of the gas, w e calculate the 
entropy accretion rate as (lGhakrabartilll9f^ . 


M 



(5) 


where, v — (7 -b l )/(7 — 1). In an accretion disk, M remain 
constant all throughout except at the shock transition where 
local turbulence generates entropy. Therefore, the entropy 
accretion rate in the post-shock region is higher than that 
in the pre-shock region. 

By definition, the accretion flow around the black holes 
must be transonic in nature. This is due to the fact that the 
velocity of the accreting matter at large distances from the 


































4 Ramiz Aktar, Santabrata Das, Anuj Nandi 


black hole is negligibly small and therefore, the flow is sub¬ 
sonic. However, flow crosses the event horizon with velocity 
equivalent to the speed of light implying that the flow is su¬ 
personic close to the black hole. This clearly indicates that 
the accretion flow necessarily changes its sonic state during 
its journey from the outer edge of the disc to the horizon. In 
order to calculate the location where the state of this sonic 
transition occurs, we derive the sonic point conditions using 
Eqs. (1-5) and obtain as. 


where. 


and 


dx^ D' 


N = 


7 -I- 1 


3 

X 


dln<I>^ 

dx 


dx 


D = V — 


2a^ 

(7 -I- l)u’ 


( 6 ) 

( 7 ) 

(8) 


where, the subscript ‘e’ signifies the quantity measured at 
the disc equatorial plane. Since the accretion flow is smooth 
everywhere, the radial velocity gradient must be finite ev¬ 
erywhere too. According to Eq. (6), if the denominator D 
vanishes at any radial distance, the numerator N also van¬ 
ishes there. Such a location, where D = N = 0, is called as 
sonic point. Setting D — 0 and N = 0 independently, we 
hnd the sonic point conditions as. 


and 


A" = 0 ^ = (7 + 1) 



1 2 

Lc — — — 

Clc 

V 7+1’ 

fd^A 

3 dlnil?^ 

\dx J ^ 

X dx 


( 9 ) 


( 10 ) 


Here, the subscript ‘c’ denotes the quantities evaluated at 
the sonic points. We use sonic point conditions, namely, Eqs. 
(9-10), in Eq. (1) to obtain the location of sonic point for a 
given set of {£, A) of the flow. For a physically acceptable 
transonic so lution, flow must contain at least one saddle type 
sonic point llPasI 1 2007ll and references therein). Depending 
on the values of £ and A, flow may possess multiple sonic 
points as well whic h is one of the ne cessary condition to 
form a shock wave ilChakrabartilll99^ . At the sonic point, 
Eq. (6) takes the form as dv/dr = 0/0 and therefore, we 
use I’Hospital rule to calculate the radial velocity gradient 
there. Once the flow variables at the sonic point are known, 
we integrate Eq. (6) starting from the sonic point inward up 
to the black hole horizon and outward to a large distance 
to obtain the full set of global accretion solution which may 
own standing shock waves. 


2.2 Governing equations for Outflow 

In this work, we consider the outflow to be emerged out 
from the accretion disc along the rotation axis of the black 
hole with the same energy and angular momentum as the 
accretion flow since we neglect the dissipative processes (i.e., 
viscosity, cooling etc.). Similar to accretion flow, we consider 
the outflow to obey the polytropic equation of state as Pj = 


Kjp'f, where, the suffix ‘j’ denotes the outflow variables. The 
equations of motion that describe the outflow dynamics are 
given by, 

(i) the energy conservation equation of outflow: 

+ TTT + ^’ 

where £j {= £) represent the specific energy of the outflow, 
Vj is the outflow velocity and aj is the sound speed of the 
outflow, respectively. 

(ii) Mass conservation equation of outflow: 


Mout = PjVjA, (12) 

where, Al is a geometrical quantity representing the total 
area function of the outflow. T o obtain A, we consi der the 
outflow geometry as described in iMolteni et al\ (Il996ll where 
the outflowing matter tends to come out through the two 
surfaces, namely the centrifugal barrier (CB) and the funnel 
wall (FW). We calculate the centrifugal barrier by identi¬ 
fying the pressure maxima surface as {d^/dx)rQg = 0 and 
the funnel wall by dehning the null effective potential as 
~ 0 dMolteni et al.|[l996l l. As the effective potential 
of our interest is complex in nature, it is unattainable to pre¬ 
vail an analytical expression of the CB surface and the FW 
and therefore, we compute them numerically. In Figure 1, 
we illustrate the outflow geometry around the rotating black 
holes for a wide range of at marked in the figure and compare 
them with the same o btained using the pseudo- Newtonian 
potential introduced bv IPaczvhski fc Wiital lll98(t) appropri¬ 
ate for a stationary black hole. The regions bounded by the 
dashed and dotted curves are for ak = —0.998 and ak = 0.8, 
respectively and the solid curves depict the result of station¬ 
ary black hole. Note that the outflow geometry for station¬ 
ary as well as for rotating black holes is indistinguishable 
outside the range of few tens of Schwarzschild radius. This 
is possibly due to the fact that the effect of black hole spin 
on the space-time geometry cursorily diminishes with the in¬ 
creasing distances. Keeping this in mind, we therefore adopt 
the outflow geometry of the stationary black hole for our 
present study in order to avoid the rigorous numerical cal¬ 
culations. This enables us to obtain an analytical expression 
of area function A and its higher order derivatives, which 
is r equired in the sonic point an alysis of outflowing mat¬ 
ter dDas fc ChattoDadhvavl [ 2 OO 8 II . A comprehensive study 
of sonic point analysis for outflows including the area func- 
tion A and its derivat i ves h ave already been presented in 
iDas fc Chattopadhyavl (l2008l f and therefore, we avoid repe¬ 
tition here. 


3 SOLUTION METHODOLOGY 

In this study, we consider the post-shock matter (PSC) as 
the precursor of the Jet base and thus, we focus only to those 
accretion solutions that possess standing shock waves. For 
shock, the accretion flow variables experience discontinuous 
transition characterize d by the Rankine-Hugoni ot (hereafter 
R-H) shock conditions (iLandau fc Lifshitzlll959ll . These con¬ 
ditions include the conservations of mass flux, energy flux 
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Figure 1. Comparison of Jet geometry for angular momentum 
A = 3.3. Funnel wall (FW) and centrifugal barrier (CB) are 
marked in the figure. Here, OC = rj is the spherical radius repre¬ 
senting the streamline of the outflow. Dashed and dotted are ob¬ 
tained for afi = —0.998 and 0.8 while solid cur ve denotes the result 
obtai ned using pseudo-Newtonian potential dPaczvnski &: Wiital 
Il980h . Jet geometries for stationary as well as rotating black holes 
are indistinguishable beyond few tens of Schwarzschild radius. See 
text for details. 


and momentum flux across the shock, respectively. In pres¬ 
ence of mass loss, a part of the inflowing matter is effectively 
emerged out as outflow from the post-shock region (PSC), 
while the remaining matter is advected in to the black hole 
straight away. Therefore, in the present scenario, the shock 
conditions are given by, 

(i) Conservation of mass flux: 

M+= M-- Maut = M-{1-Rrh). (13a) 

The quantities having subscripts and ‘-I-’ are referred 
to the values before and after the shock which we follow 
throughout the paper unless otherwise stated. The pre-shock 
and post-shock accretion rates are denoted by M_ and M +, 
respectively and Mout is the mass flux for the outflowing 
matter. Following this, the mass outflow rate is computed 
as Rrh = Mout/M-. The next condition is 

(ii) Conservation of energy flux: 

£+=£-, (136) 

and finally, we have 

(iii) Conservation of momentum flux: 

W+ + 'E+vl = W- + T.-v'^, (13c) 

where, W an d S represent the vertically integrated pressure 
and density llPas et a6ll2001bll . 

We rewrite Eq. (13b) and Eq. (13c) in terms of Mach 
number {M = vja) and are given by, 


and 


2a+ 

(37 - l)M+ 




l-i?ri 


2a_ 

(37 - 1)M_ 


-\- NT— a_ 


_ (13e) 

where we utilize Eq. (13a) and a = ^(37 — 1)IE/(2E) 
dDas et al. I l2001bll . 

Using Eqs. (13d-e), we obtain a shock invariant quantity 
(Cs) in terms of Mach number M as, 


Cs = 


[ iv ^ + (37-1)M+]'(1-7?. 


[i^ + (37-1)M-]' 


[2 -b (7 - 1 )M 2 ] [2 + (7 - l)Ml] 

where, M- and M+ stand for Mach number just before and 
after the shock in the pre-shock and post-shock flow, respec¬ 
tively. 

Since the shock conditions are coupled with inflow and 
outflow variables, the accretion and Jet equations are solved 
simultaneously. We obtain the Jet velocity gradient from its 
governing equations and calculate the s onic point propertie s 
following the ‘critical point’ analysis (IChakrabartil Il989l l. 
While doing this, we use the inflow parameters, namely, 
energy £ and angular momentum A of the inflow. This is 
because the outflow is considered to be originated from the 
the post-shock region (PSC). This also ensures that the Jets 
are launched with the same density as the post-shock flow, 
namely pj = p+. Further, we integrate the Jet equations to 
calculate the outflow variables at the Jet base starting from 
the Jet sonic point and using Eqs. (3), (4) and (12) we com¬ 
pute the mass outflow rate in terms of the inflow-outflow 
properties at the shock and is given by. 


^ _ Movf _ PjVj{Xs)A{Xs) _ RVj(Xs)A{Xs) 

^ ^ ATVP-V-Xsh- ^ A7v^xl^^K~"^^a+v- ’ 

(14) 

where, R — E+/E_(= p+6,+ /p_/i_ = pja+/p-O-), is the 
compression ratio and Vj{xs) and 7l(a;s) are the Jet velocity 
and the Jet area function at the shock, respectively. We 
utilize an iterative method to calculate Rm self-consistently 
which is given as follows: 

We begin with Rrh = 0. Using the shock invariant quantity, 
we calculate the virtual shock location for a given set of 
{£, X) with the consideration that the entropy of the inflow¬ 
ing post-shock matter (A4 + ) and outflowing matter (Mout) 
are larger than the entropy of the pre-shock matter (Ad-). 
This eventually reflects the fact of second law of thermody¬ 
namics as the shocked solutions ascertain the preferred mode 
of accretion. We use the same set of inflow parameters for Jet 
equations and calculate Rrh which we employ further in the 
shock invariant equation to obtain a new shock location. We 
continue this iteration process until the solution converges 
to the actual shock location and accordingly we compute the 
corresponding Rrh- In the following Sections, we investigate 
the properties of Xg and Rhn in terms of the inflow param¬ 
eters {£, A) for various values of black hole spin (at) and 
present the results which are followed by the discussion on 
astrophysical applications. 


4 RESULTS AND DISCUSSIONS 


1 ,,r2 2 , 

_M+a+ + ^ 


= (13d) 

2 7 “ t 


In an advective accretion process around black holes, the ac¬ 
creting matter suffers discontinuous shock transitions due to 
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the centrifugal barrier. This causes the post-shock flow to be¬ 
come hot and puffed up which eventually behaves like a effec¬ 
tive boundary layer around the black hole (i.e., PSC). Dur¬ 
ing accretion, a part of the inflowing matter is deflected ver¬ 
tically by this boundary layer to produce thermally driven 
outflows. To illustrate the complete picture of the accre¬ 
tion and ejection mechanism, we present the global inflow- 
outflow solutions for various at in Fig. 2. In panel (a), (b) 
and (c), we show the variation of Mach number {M = v/a) 
of the inflow with the radial distance (x) for ak = 0.8, 0.0, 
and -0.998, respectively. Subsonic matter at large distance 
gradually gains its radial velocity as it proceeds towards the 
black hole due to the fatal attraction of gravity and crosses 
the outer sonic point (xout) to become supersonic as shown 
by the arrows. As the flow continues its journey further to¬ 
wards the black hole, it experiences discontinuous transi¬ 
tion in flow variables when the R-H conditions are favorable 
which is indicated by the dashed vertical arrow. After this 
transition, flow momentarily slows down and gradually picks 
up its velocity. Eventually, flow becomes supersonic again 
after crossing the inner sonic point (xin) and finally enters 
in to the black hole. As the outflow is emerged out from 
the effective boundary layer of the black hole, the density 
of the post-shock flow is decreased causing the reduction 
of pressure at PSC as well. In order to maintain the pres¬ 
sure balance across the standing shock in presence of mass 
loss, shock itself has to move inward which is indicated by 
the solid vertical arrow. In the panel (d), we plot the Mach 
number {MR variation of the outflowing matter with its ra¬ 
dial coordinate {xR corresponding to panel (a), (b) and (c). 
The black solid circles represent the sonic locations and the 
arrows show the direction of motion of the outflowing mat¬ 
ter. The inflow parameters for panel (a) is (T, A) =(0.007, 
2.65), for (b) is (0.0025, 3.45) and for (c) is (0.0015, 4.0), 
respectively. 

In Fig. 3, we compare the location of the shock transi¬ 
tions as function of energy S for a set of angular momentum 
A. In the upper panel, we choose ak = 0.6 and find that 
the inflow-outflow solution possesses shock wave for a wide 
range of S and A. As anticipated in Fig. 2, the shock location 
moves towards to the black hole horizon when Rth R 0 just 
to maintain the pressure balance across the shock. We again 
observe that the shock location reaches to its lowest value 
around x™‘" ~ 8rg at energies higher compared to the case of 
Rrh ~ 0 and the limiting value of this energy increases with 
the decrease of angular momentum of the flow. However, 
the maximum value of energy that allows shock transition 
in presence of mass loss is identical with no mass loss case. 
This provides a clear hint that the range of inflow parame¬ 
ters for outflows, namely energy £ and angular momentum 
A are reduced from their lower ends. In the lower panel, we 
consider ak — —0.6 and obtained the similar results which 
differ only quantitatively. Overall, it is clear that the possi¬ 
bility of shock formation is affected substantially due to the 
presence of outflow. 

Before we proceed further, we now investigate the effect 
of the black hole rotation on the generation of mass outflow 
rate. In Fig. 4, we present the variation of shock location 
(xs) and outflow rate (i?m) as function of a*,. Here, we fix 
£ = 0.002 and vary A from 2.8 to 3.92 from the right most 
curve to the left most with an interval AA = 0.16. In the 
upper panel, the solid curves represent the shock locations 
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Figure 2. Variation of the inflowing Mach number M (= v/a) 
with radial distances (x) for (a) ak = 0.8, (b) ak = 0 and (c) ak = 
—0.998, respectively. The corresponding outflow Mach number 
Mj{= Vj/aR variation is shown in panel (d). See text for details. 




Figure 3. Variation of shock location as function of energy £. 
The solid curves represent results without mass loss and the dot¬ 
ted curves are with mass loss. In the upper panel, we choose 
ak = 0.6 and curves are for A = 2.78 (right), 2.80 (middle) and 
2.82 (left), respectively. In the lower panel, we consider ak = —0.6 
and curves are for A = 3.68 (right), 3.70 (middle) and 3.72 (left), 
respectively. 
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Figure 4. Variation of (a) shock location {xs) and (b) outflow 
rate (Rrh) 3-s function of black hole rotation parameter aj. for A 
= 2.8 to 3.92 (right to left), where AA = 0.16. Here, we fixed 
the flow energy as f = 0.002. In the upper panel, dashed curves 
denotes the variation of shock location in presence of outflow. 



Figure 5. Variation of (a) shock location (xs), (b) outflow rate 
{Rrh) and (c) compression ratio {R) as function of energy € for 
aj, = 0.6,0.65,0.7,0.75 and 0.8 (right to left). Here, we fix the 
flow angular momentum A = 2.8. 


in absence of mass loss. However, in presence of mass loss, 
accretion flow adjusts the location of the shock transition 
closer to the black hole horizon in a way that the pressure 
balance condition across the shock is maintained which is 
depicted by the dashed curves. On the contrary, as at is 
increased for a set of £ and A, the shock locations recede 
away from the black hole horizon. This indicates that the size 
of the post-shock region (PSC) is enhanced with the increase 
of ak and the inflowing matter is eventually intercepted by 
the large effective area of the post-shock flow that produces 
more outflow rate. In the lower panel, we present the feature 
of outflow rate variation with a*,. For a set of £ and A, the 
outflow rate Rm shows non-linear correlation for prograde 
as well as retrograde flows. In addition, we observe that for 
a given a*,, Rm is higher for increasing A and with this, we 
infer that large outflow rate is associated with the higher A 
and lower ak when £ is fixed. Overall, we find that outflow 
rate reaches close to 16% for these parameters as depicted 
in the figure. 

Next, we study the characteristics of the post-shock 
quantities in terms of the inflow variables and present them 
in Fig. 5. Here, we choose the angular momentum as X — 2.8 
and vary energy £ and black hole rotation parameter ak¬ 
in the upper panel, we plot the variation of shock location 
(Xa) with £ for ak varied from 0.6 (right) to 0.8 (left) with 
Aafe = 0.05. The corresponding variation of outflow rate Rm 
is shown in the middle panel. Due to shock transition, the 
post-shock flow is compressed. The measure of this compres¬ 
sion is quantified as the ratio of the post-shock density to 
the pre-shock density and it is termed as shock compression 
ratio R. In the lower panel, we present the variation of R 
with £. When shock forms closer to the black hole horizon. 


the amount of gravitational potential energy release is higher 
ensuing the formation of strong shock. As the energy is in¬ 
creased, the shock location moves outward that increases the 
outflow rate and weakens the shock compression ratio. We 
observe similar behaviour for all ak- However, for a given A 
and ak, there is a range of energy beyond which mass out¬ 
flow ceases to exist. This provides an indication of finding 
a parameter space for shock in the fundamental plane of 
energy and angular momentum of the inflow. We perform 
the similar study for retrograde flow as well and find the 
identical trend as depicted in Fig. 6. 

As we pointed out earlier that the transonic accretion 
solutions including R-H shock waves are not isolated solu¬ 
tions. In fact, these solutions do exist for a wide range of 
parameters, namely, energy {£), angular momentum (A), 
and black hole rotation parameter {ak). Towards this, we 
identify the shock induced global accretion solutions using 
the parameter space spanned by the energy (£) and angular 
momentum (A) of the flow and classify them in terms of ak- 
In Fig. 7, we separate the parameter spaces with the solid 
boundaries that indicate the regions for global shock accre¬ 
tion solutions in absence of mass loss. The corresponding 
values of ak are marked in the figure. We further investi¬ 
gate the parameter spaces that cater mass loss and indicate 
it with the dashed boundaries. Due to mass loss, the post¬ 
shock region (PSC) shrinks as seen in Fig. 2, and therefore, 
the associated parameter space is reduced compared to the 
results having no outflows, particularly towards the lower 
energy and lower angular momentum sides (see Fig. 3). For 
flows with input parameters chosen from these part of the 
parameter space, shock forms v ery close to the black hole 
horizon when outflow is ignored llPas et al ]|2001bl) and they 
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Figure 6 . Same as Fig. 5. Here, = —0.998, —0.8, —0.6, —0.4, 
and —0.2 (right to left) and flow angular momentum is A = 3.9. 



2 2.4 2.8 3.2 3.6 

X 


Figure 8 . Effective regions of parameter space separated in terms 
of shock waves transition for various values of 7 marked in the 
figure. Here, we choose = 0.5. The solid boundary denotes the 
presence of shock waves in absence of mass loss and the dashed 
boundary represents the results including mass loss. 



Figure 7. Energy angular momentum parameter space. Solid 
curve (black): shock parameter space without outflow. Dashed 
curve (red): shock parameter space with outflow. Various values 
of are marked in the figure. 


cease to exist when the outflow is allowed to emerge out 
from the post-shock disc. 

So far, we have presented the shock induced complete 
global inflow-outflow solutions and its properties for a given 
value of 7 = 4/3. In reality, however, the theoretical limit 
of 7 lies in the range between 4/3 to 5/3 depending on 
the ratio between the thermal energy and the rest energy 


of the flow ([Frank et a/.ll20Q3 ). To infer this, we consider 
rotating flows that are characterized as thermally ultra- 
relativistic (7 ~ 4/3), thermally trans-relativistic (7 ~ 1.4) 
and thermally semi -non-relativistic (7 ~ 1-5), respectively 
dKumar et al ][ 2 oH) and obtain the parameter spaces for 
Ofe = 0.5 similar to Fig. 7. This is shown in Fig. 8 . Here, 
we observe that in all three cases shocks exist for significant 
range of inflow parameters. However, the range of parame¬ 
ters for shock reduces considerably as the flow changes its 
ultra-relativistic character towards the non-relativistic limit. 
In addition, we find that all the parameter spaces further 
shrink as in Fig. 7, when a part of the inflowing matter is 
deflected in the form of mass loss from the inner part of the 
disc (PSC). 

According to the formalism adopted in this work, the 
mass outflow rate Rrh is computed self-consistently in terms 
of the inflow parameters around a black hole of rotation pa¬ 
rameter Ofe. This allows us to estimate the maximum R^°‘^ 
for a given value of 7 . While doing this, we identify a par¬ 
ticular set of energy {£) and angular momentum (A) from 
their full range that provides the highest value of Rm- The 
importance of this investigation is associated with the study 
of the maximum Jet kinetic power corresponding to 
which we discuss in the next section. Following Fig. 8 , we 
calculate the variation of as function of au for various 

values of 7 which is presented in Fig. 9. The results dis¬ 
played in upper, middle and bottom panels are for 7 = 4/3, 
1.4 and 1.5, respectively. We And that gradually in¬ 

crease with the increase of a* for all cases. We also observe 
that the accretion flows having 7 = 4/3 have the poten¬ 
tial to produce more outflows compared to the flows with 
higher values of 7 . This is perhaps inevitable due to the fact 
that the outflows under consideration are thermally driven 
and therefore, thermally ultra-relativistic flows exhibit max- 
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Figure 9. Variation of maximum outflow rates with the 

black hole rotation parameter a^. Upper panel is for 7 = 4/3 (a), 
middle panel is for 7 = 1.4 (b) and bottom panel is for 7 = 1.5 
(c), respectively. See text for details. 



tion states with time and the characteristic features of the 
emergent radiations and the Jet properties are also var¬ 
ied accordingly. During these dynamical change of states, 
it is unlikely that the dissipative properties (i.e., viscos¬ 
ity, cooling processes etc.) of the accretion flow will re¬ 
main constant all throughout, instead, they seem to adjust 
in a such way that suitably represents the dynamical vari¬ 
ation of PSC. In LHS of GBHs, the typical geometry of 
PSC is quite large which reduces considerably in HIMS as 
the dynamical shock moves towards the black hole during 
the rising phase of the outburst. Such a trend has been re¬ 
ported while modeling the e volution of QPO frequen c ies for 
several black hole sources (|Oialmibaj1J__ej_aLl 200^ 200^ 
Nandie|_a^ 2oS_ Debnath et al\ 20131 : Radhika fc Nandi 


2014 liver 6101 . 11201 !^ . Interestingly, it is also observed that 


most of the GBH sources show per sistent/steady radio emis - 
sion during these accretion states dPender et oZ. [2004 . l2009l i. 
Now, we attempt to address the evolution of accretion states 
from LHS to HIMS and its association with mass loss using 
our present formalism as case by case which possibly acts 
as a local model at individual time frame. Toward this, a 
set of energy and angular momentum {£, A) of the inflowing 
matter is identified which would represent their local values 
for a dissipative accretion flow and obtain the shock location 
Xs and outflow rate Rm for a given a^. Since the post-shock 
geometry (PSC) is associated with Xs, as anticipated ear¬ 
lier, we chose various such sets of {£, A) in a way that Xa 
moves towards the black hole and show the variation of 


with Xs in Fig. 10 for a typical value of a*, = 0.5. Note that 
Rrh is increased with the decrease of Xs resulting more Jet 
kinetic power that represents the common characteristic of 
radio emissions observed in outbursting sources. 


5 ASTROPHYSICAL APPLICATION 

Here, we attempt to estimate the mass outflow rate and its 
associated Jet kinetic power for various astrophysical black 
hole sources (GBHs and AGNs) using our present formal¬ 
ism. Since our work deals with the steady outflows, we focus 
only to those sources particularly to their accretion states 
where persistent/steady Jets are observed. These Jets are 
essentially compact (i.e., optically thick) in nature and are 
not isolated completely from the core of the central engine 
llMirabel fc Rodrigue Jl 19^ . Il998l : iGorbel et oZ.ll2003 . 12001 I : 
iFender et aZ.ll200lh. 


Figure 10. Variation of outflow rate (Rrh) with shock location 
{xs) obtained for different sets of energy {£) and angular momen¬ 
tum (A). Here, aj, = 0.5 is used for representation. See text for 
details. 

imum outflows. Following this, in the next Section, we con¬ 
sider 7 = 4/3 while computing the Jet kinetic power for 
GBHs and AGNs, until otherwise stated. 

Until now, we discuss the various properties of the 
accretion-ejection solutions in terms of the inflow parame¬ 
ters. However, all these studies are based on the assump¬ 
tion of stationary state that does not represent the dy¬ 
namical behaviour of accretion flows around the outburst¬ 
ing BH sources. Usually, these sources change their accre- 


5.1 Jet kinetic power predicted from Accretion 
States 

We consider five black hole sources namely, XTE 
J1859-f226, GRO J1655-40, GX 339-4, H 1743-322 and GRS 
1915-1-105, respectively. For these sources, mass (Mbh), dis¬ 
tance (d) and spin (at) are constrained within the accu¬ 
racy limit and are given in Table 1. In order to estimate 
the Jet kinetic power, we calculate the X-ray flux for these 
sources in the low-hard state (LHS) and hard-intermediate 
state {HIMS) using RXTE^ satellite data. 


In general, GBH sources undergo ou t bursts 


llHoman & Bellonil l2005l: Remillard & McGlintockI 

2006|: 

Debnath et aZ.ll2008l:lNandi et oZ.ll2012l:lDebnath et al\ 

2 OI 3 I: 

Radhika & Nandil 2014: Tver et aZ.ll2015f) and show activitv 
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of radio emissions coupled with the accreti on states of their 
evolu t ion during the ou t burst phases (iBrocksopp et gZ] 
2OO2I: iFender et^ I2OO4I. [20nil: ICadolle Bel et al\ I2OIII: 
Miller-Jones et al\ I2OI2I : iRadhika fc Nandil 2014 1. It has 
also been observed that most of the outbursting sources 
show persistent/steady Jet activity during LHS and HIMS 
while transient ‘relativistic’ ejections are observed during 
the tra nsition from hard-intermediate to soft - intermediate 
state l| Brocksopp et al. I I2OO2I: iFender et al. I I20o4 I2009I: 
IRadhika fc Nan^ 2014 1. 

In order to estimate the ‘model predicted’ Jet power 
based on accretion states, we estimated X-ray flux from a 
single observation of each states (i.e., LHS and HIMS) for 
all the sources from different outburst phases (as mentioned 
in the Table-1). We use archival data obtained from the 
HE AS ARC database of the RXTE satellite for the estima¬ 
tion of X-ray flux for all the five sources. We extracted and 
analyzed background-subtracted PCA (3 - 30 keV) spec¬ 
tral data using PCU2 detector (i.e., well calibrated detector) 
for our specific purpose. The standard FTOOLS package of 
HEASOFT v6.15.1 was used for spectra l data reduction (see 
iNandi et al\ 1I2OI2II : IRadhika fc Nandi ll2014l l for details). 
For spectral analysis and modeling, we used the packages 
of XSPEC V12.8.1. 

We have done spectral modeling for each observations 
from each spectral states (i.e., LHS and HIMS) for all the five 
sources using PCA spectral data in the energy range of 3 - 30 
keV. We model the energy spectrum using the phenomeno¬ 
logical accretion disk model i.e., consisting of a diskbb and 
a power-law component. In this modeling, the diskbb and 
power-law components provide the contribution from the 
accretion disk and ‘hot’ Compton corona (i.e., PSC) and 
thereby one can estimate the total unabsorbed X-ray flux 
{Fx, in units of ergs cm~^ s~^) emitted from the accretion 
disk around the GBH sources. 

Once we obtain the unabsorbed X-ray flux (Fx) of a 
source, we calculate the X-ray luminosity (Ta,) of the source 
employing the relation L^ = Andi^F^, where d is the distance 
of the source. Assuming the maximum radiative efficiency of 
the infalling matter around rotating black hole is 7; ~ 0.3, 
we calculate the accretion rate of the black hole as, 

= 2.99 X 10-« 

Here, Mbh denotes the mass of the black hole and MEdd 
represents the Eddington accretion rate. Next, we calculate 
the mass outflow rate using our theoretical estimate (Eq. 
14) which has its maximum limit (J?™’^”') for a particular set 
of {ak,£,X) (see Fig. 9). With this, we find the maximum 
mass outflow rate as Mout = BRNMace which successively 
allows us to compute the maximum Jet kinetic power as, 

L max 7 -»max ..nV .,2 —1 /-i r\ 

Jet = Rra X Macc X c ergs s . (15) 

As pointed out earlier that steady and pe rsistent radio 
emissions a re observed during LHS and HIMS ([Fender et al\ 
I2OO4I . I^Ogh . In addition, there are indications that the ra¬ 
dio emission increases while the state transition from LHS 
to HIMS takes place. Moreover, the radio emission becomes 
non-steady during the HIMS itself just prior to the ejection 
and it is quite known that the relativistic ejecti ons take place 
durin g the transition from HIMS to SIMS (IFender et al\ 
I2OO4I . ]2009h . All these observations perhaps indicate that 


the persistent Jet activity would be maximum during HIMS. 
Therefore, in order to study the outflow properties in the 
HIMS, we consider the maximum outflow rate obtained from 
our model calculation (see Fig. 9) and use it to estimate the 
Jet kinetic power. On the other hand. Jet kinetic power in 
the LHS is expected to be lower than the HIMS although its 
quantitative estimate is unclear. Therefore, we use outflow 
rate to be around ~ 10% as a representative value while cal¬ 
culating the Jet kinetic power for LHS. Below, we mention 
the fundamental properties of each sources and present the 
details for the estimation of Jet kinetic power corresponding 
to individual states of a particular observation. 


XT E J1859 -h 226: _ 

iFilippenko fc ChornockI (|200ll 'l first presented the dy¬ 
namical estimate of m ass of the source to b e around 
7.4 ± I.IMq. Recently, IRadhika fc Nandil (l2014l l claimed 
that the mass of XTE J1859 + 226 is perhaps in between 
6.58M(Vi — 8.84M(Vi which is sim ilar to the prediction of 
IShaposhnikov fc TitarchukI (l2009l) although the lower mas s 
limit is estimated as 5.4M0 bv ICorral-Santana et al\ (l201llf . 
However, we consider the typical mass of the source as 
7Mr^. The distance of this s o urce is around d ^ llkpc 
([Filippenko fc ChornockI 1200 ill . ISteiner e t al\ (|2 013tl mea- 
sured the spin as a^ ~ 0.4, however, Motta et al\ l2014al l 
recently reported the spin of the source is au ~ 0.34. 
Since the spin predictions are quite close, we use a^ ~ 0.4 
for this analysis. We estimate the fluxes Fa, (see Table 
1) of LHS and HI MS of the 2009 outburst of the source 
(IRadhika fc Nan^l2014ll . The corresponding disc luminosi¬ 
ties are calculated as L^^^^ = 8.26 x 10®^ ergs s“^ and 
= 1-85 X 10^® ergs s“^, respectively. Now, it is 
reasonable to assume the accretion efficiency for rotating 
black hole as r; = 0.3 which corresponds to the accretion 
rate of the inflowing matter as Ma^^ = Q.'XtAMEdd in 
LHS and = 0.680MBdd in HIMS. For LHS, we 

use Rm = 9.83% following our theoretical estimate where 
Xa = 64.6rg for Ok = 0.4, £ = 0.00198 and A = 3.18. Incor¬ 
porating these inputs in Eq. (15), we obtain the Jet kinetic 
power as Lj^t^ = 2.52 x 10®^ ergs s“^. The maximum mass 
outflow rate for HIMS corresponding to at = 0.4 is obtain 
from Fig. 9 as R)^^^ = 17.5% for £ = 0.00547, and A = 3.1, 
where the shock transition occur at 21.9rg. Using these val¬ 
ues in Eq. (15), we obtain the maximum Jet kinetic power 
as Lj^A^^ = 1.08 X 10®® ergs s“^ which we regard to be 
associated with the Hard-Intermediate state of this source. 


GRO J1655-40: 

The mass of the source GRO J1655-40 is reported 
by Greene et al. ( 200 il l and is given by 6.3Mq. Recently, 
iMotta et al. ( 2014bl 'l estimated the object mass as 5.31 ± 
0.07 Mq. This source is located at around d ~ 3.2 kpc 
and t he spin of this source is estimated by IShafee et al\ 
(l2006tl using RXTE^ and ASCA data through the model¬ 
ing of the thermal spectral continuum and obtained in the 
range a*, ~ 0.65 — 0.75 (and reference therein). Fitting the 
strong reflection features of iron line in XMM-Newton data, 
iReis et ^ (120091 ') determines the lower limit of the spin of 
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Table 1. Accretion state dependent Jet kinetic power. 


Objects 

Mbh 

d 


Observation (3 — 30 keV) 

^acc 

Rrh 

Rjet 

rObs # 

^Jet 


(Mq) 

(kpc) 


States 

Fx (ergs cm ^ s ^) 

(MEdd)* 

(%) 

(ergs s~^) 

(ergs s~^) 

XTE J1859+226 

7 

11 

0.4 

LHS 

5.71 x 10“® 

0.304 

9.83 

2.52 X 10®’^ 


(2009 Outburst) 




HIMS 

12.79 x 10“^ 

0.680 

17.5 

1.08 X 10®® 

1.82 X 10®® (1) 

GRO J1655-40 

6.3 

3.2 

0.7 

LHS 

3.19 x lO-'^ 

0.016 

9.98 

1.18 X 10®® 


(2005 Outburst) 




HIMS 

8.02 x 10“® 

0.040 

17.68 

5.79 X 10®® 

3.01 X 10®® (2) 

GX 339-4 

7.5 

15 

0.4 

LHS 

15.71 x lO-** 

1.450 

9.98 

1.29 X 10®® 


(2002 Outburst) 




HIMS** 

12.12 x 10“® 

1.118 

17.5 

1.90 X 10®® 

2.92 X 10®® (1) 

H 1743-322 
(2009 Outburst) 

8 

8.5 

0.2 

HIMS 

4.03 x 10“® 

0.112 

17.35 

2.01 X 10®’’ 

1.08 X 10®® (3) 

GRS 1915-1-105 
(1997 Observation) 

12.4 

8.6 

> 0.98t 

LHS 

20.33 X 10"® 

0.373 

9.98 

5.99 X 10®’' 

8.06 X 10®’’ (1) 


* MEdd = 1-44 X 10^’^ gm s 1 

** Total flux in HIMS (3 — 30 keV) is smaller than LHS as the contribution of the hard X-ray flux (> 10 keV) in HIMS is less. 

# LObs = ^Moutc2 is used for the source H 1743 — 322, where Mout is the outflow rate (see reference). For other sources, = Lr X LEdd 

is used, where Lr is the observed Jet power (in Edd) (see references) and LEdd = 1-3 x 10®® ( s“®. 

References: fli lFender et all ||2004| ~) (2) iMigliari et oD ll2007l ) (3i iMiller et al\ ll2012li 


this object at = O.Q. lMotta et al\ ll2014hJ ) calculated the spin 
of the object using X-ray timing method and found as = 
0.290 ± 0.003 which is an inconsistent estimate compared 
to iron line or continuum methods. In our present analysis, 
however, we consider a*, = 0.7. As before, we analyzed the 
2005 outburst of GRO J1655-40 to calculate the X-ray fluxes 
(Fx) for LHS and HIMS (see Table 1) in the energy range 
3 — 30 keV and obtain the corresponding accretion rates 
= 0.016M_Bdd in LHS and = 0.040M_Bdd 

in HIMS. We use Rrh = 9.98% in case of LHS which is 
obtained for = 0.7, S = 0.00258 and A = 2.845 where 
Xs = 49.98rg. From these values, we estimate Jet kinetic 
power in LHS as LjS® = 1-18 x 10®® ergs s“®. For HIMS, 
we estimate of the maximum mass outflow rate for = 0.7 
(see Fig. 9) as = 17.68% where £ — 0.0073 and 

A = 2.715 with Xs = 16.73rg. These inputs provide the max¬ 
imum Jet kinetic power as = 5.79 x 10®® ergs s“®. 

GX 339-4: 

_ T he ma ss of the object is estimated as 7. 5 ± 0.8 Miv, 

[Chen (l201ll ) and distance d ~ 15 kpc by iHvnes et al\ 
1120041) . The issue of spin measurement of this objects is 
not settled yet as there are conflicting measurements. An¬ 
alyzing the XMM-Newto n dat a set for bro a d iron line de¬ 
tection, iReis et ^ (l2008l ) and iMiller et al\ ll2008l) claimed 
flfc = 0.935 while the disc continuum f itting prefers the lower 
spin having upper limit of ak < 0.9 dKolehmainen fc Pond 
2010|). Howeve r , ana lyzing the wide-band Suzaku speatra, 
Yamada et al\ (l2009l) pointed out < 0.4. Being aware 
of these, we consider a conservative estimate of spin as 


ak = 0.4 in our calculation. This source has undergone out¬ 
burst phases several times during RXTE era. For this anal¬ 
ysis, we consider the 2002 outburst spectral data (in the 
energy band of 3 — 30 keV) for LHS and HIMS. During May 
2 of 2002 outburst, the object was in LHS emitting X-ray 
flux of = 15.71 X 10“® ergs cm“^ s“^ (see Table 1) 

and the disc luminosity of L^Hf = 4.22 x 10®® ergs s“®. 
Using Tj = 0.3, the disc rate is calculated as = 

1.450M_Bdd- Estimating Rdi = 9.98% from our model us¬ 
ing £ = 0.00198 and A = 3.18 with Xs = 64.63rg, we cal¬ 
culate the Jet kinetic power — 1.29 x 10®® ergs s“®. 

On May 12 of the 2002 outburst, the source was in HIMS 
and the radiated X-ray flux was calculated as F^^^^ = 
12.12 X 10“® ergs cm~^ s“® (see Table 1) and the disc lu¬ 
minosity of = 3.26 X 10®® ergs s“®. Using rj = 0.3, 

the disc rate is calculated as = 1.118f\j^Bdd. Esti¬ 
mating = 17.5% from our model using £ = 0.00547 

and A = 3.1 with Xs = 21.9rg, we calculate the maximum 
Jet kinetic power = 1.90 x 10®® ergs s“®. 

H 1743-322: 

IMiller et al\ ||2012| ) and ISteiner et al\ (l2012l) reported 
the mass, distance and spin of H 1743-322 as ~ 8Mq, 
d ~ 8.5 kpc and at = 0.2, respectively. Using the RXTE 
observation of the 2009 outburst of the source, the disc lu¬ 
minosity for this source is computed for the energy range 
of 3 — 30 keV corresponding to the measured X-ray flux of 
pHiMS ^ 4 Q3 ^ ;^q-9 ^^-2 g-i Table-1) 

llMiller-Jones ef^l2012l ) as 4.63 x 10®^ ergs s“®. Assum¬ 
ing T] = 0.3, we obtain the accretion rate as = 
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Q.112MEdd- Accretion rate in LHS is not known due to lack 
of ob servation during the 2009 outburst (iMiller-Jones et ~al\ 
l2012h . For ufc = 0.2, E — 0.00485 and A = 3.176, we obtain 
the maximum mass outflow rate for HIMS as 7?™“ = 17.35% 
with Xs = 23.66rg. Employing these values, we find the max¬ 
imum Jet kinetic power as = 2.01 x 10®^ ergs s“^. 

GRS 1915-fl05; _ 

Earlier, iGreiner et al\ (1200 il l estimated the mass of 
GRS 1915-1-105 as (14 ± 4) Mq. Recently, iHurlev et gZ.I 
ll2013lf reported the new mass estimate as (12.9 ± 2.4) M(^ 
which is further revised as (12.4 ± 2) Mq by iReid ef'51 
ll2014ll . The dis tance of the source is reported as d ~ 
(9.4±0.2) kpc bv iHurlev et al\ (l2013l ) and iReid et al\ (l2014l f 
claimed the distance to be d ~ (8.6 ± 2) kpc. However, 
in this calculation we use the mass and distance of the 
source as 12.4Mq and 8.6 kpc, respectively in order to 
compute the disc luminosity. The source is extremely ro¬ 
tating as iMcClintock et oTI ll2006tl estimated the spin pa- 

rameter > _ 0.98 which is similar to the estimate of 

iBlum et al\ (l2009l ). As the source is highly variable in X- 
rays, we choose one observation when the source was in the 
hard state that is similar to LHS of other GBH sources. 
RXTE observed the source in the energy range of 3 — 30 
keV on 22*^ October, 1997 and object was in so-called x~ 
class (iBelloni et oZ.|[200ll : lNandi et oZ.ll2001cl f . For this obser¬ 
vation, the disc luminosity corresponding to the X-ray flux 
of = 20.33 X 10~® ergs cm“^ s“^ (see Table 1) is ob¬ 

tained as 1.79 X 10^® ergs s“^ which provides the disc accre¬ 
tion rate as = 0.373MEdd- According to our model, 

the mass outflow rate for LHS (or simply hard state for this 
source) is calculated as Rm = 9.98% where, Xs = 50.47rg 
with ak = 0.98, £ = 0.00276 and A = 2.679. Here, we cross 
the upper limit of ak (indicated by dagger (^) in column 4 
of Table 1) as the adopted black hole potential satisfactorily 
describes the space-time geometry for ak < 0.8. However, 
we anticipate that the obtained Rm provides qualitative 
estimate that would not differ significantly from its exact 
value. Using these values, the Jet kinetic power is found to 
be = 5.99 x 10®’^ ergs s“^. 

In this Section, we calculated the Jet kinetic 
power mostly for outbursting GBH sources (except GRS 
1915+105) for different accretion states (i.e., LHS and 
HIMS). Our findings clearly show that as the sources transit 
from LHS to HIMS, there is significant increase in the Jet 
kinetic power. The predicted Jet kinetic powers are in close 
agreement with the observed values for XTE J1859+226, 
GRO J1655-40, GX 339-4, GRS 1915+ 1 05 and H 1743-32 2 
dFender et a7|[20^ : iMigliari et oZ.II^OOTI : iMiller et a/.ll201^ ') . 


5.2 Jet kinetic power estimated from Accretion 
Rates 


We further extend our study to estimate the Jet ki¬ 
netic power for other GBH and AGN sources where we do 
not observe fast (~ day scale) state transitions similar to 
outbursting GBH sources (except XTE J1550-564). In or¬ 
der to calculate the Jet kinetic power for these sources, we 
find their mass (Mbh), accretion rate (Macc) and spin (a*,) 
from the existing literature. Meanwhile, we compute 77)0“’' 


corresponding to black hole spin ak using our theoretical 
approach. We then use the values of Macc and in Eq. 

(15) to estimate the Jet kinetic power. 

In Table 2, we display the physical parameters of the 
sources (GBHs and AGNs) under consideration along with 
the computed Jet kinetic power obtained from our analy¬ 
sis. In column 1-4, we present the list of sources, their mass 
(Mbh), accretion rate {Min) and spin {ak). In column 5-6, 
we mention the representative values of energy £, angular 
momentum A of the inflow that provide the shock location Xs 
(in column 7) and the corresponding maximum mass outflow 
rate R’^^^ (in column 8). Finally, in column 9, we present the 
maximum Jet kinetic power 71™“. The first seven sources 
are GBHs whereas the last three sources are AGNs. Here, 
we give emphasis on the maximum Rrh in order to specu¬ 
late the upper limit of Jet kinetic power. We find that the 
estimated Jet kinetic powers are in close agreement with 
the observed values at least for few so urces, namely Gyg 
X-1, XTE J1550-564, M87 and Sgr A* llFender et adlioMr 
Ide Gasperin et ^l2012l : iFalcke fc Biermann 1 19991) . For re¬ 
maining sources, we argue that the present method illus¬ 
trates the typical estimates of Jet kinetic power that possi¬ 
bly lie within the acceptable range. Further, in our analy¬ 
sis, we choose four sources having black hole spin ak > 0.8 
which are indicated by dagger (^) in column 4. As before, we 
anticipate that the obtained 71™“ for these sources provide 
qualitative estimates that would not differ significantly from 
their observed values. 


6 CONCLUDING REMARKS 

In this work, we self-consistently examine the accretion- 
ejection mechanism around the rotating black holes. We 
consider the structure of the accretion disc to be stationary, 
thin, rotating and advection dominated which contains R-H 
shock waves. These shocks are formed as a consequence of 
centrifugal barrier located near the black hole horizon. Based 
on the second law of thermodynamics, we argue that the 
shock solutions are dynamically preferred over the smooth 
solutions as they possess higher entropy (Becker & Kazanas 
2001). During accretion, a part of the super-sonic matter 
is deflected at the centrifugal barrier due to excess thermal 
pressure caused by the shock compression and eventually 
emerge out in the form of thermally driven outflows. These 
outflows are further channeled through the confined geom¬ 
etry bounded by the funnel wall and pressur e maxima sur¬ 
face along the rotation axis of the black hole llMolteni et al\ 
Il996h . 

We And that mass loss can occur for prograde as well as 
retrograde flows. When outflow is emerged out from the in¬ 
ner part of the disk (PSG), post-shock pressure is decreased 
that essentially compels the shock front to move forward 
towards the horizon in order to maintain the pressure bal¬ 
ance across the shock. Therefore, a flow originally contain¬ 
ing shock wave close to its minimum location (a;™™) for 
77m = 0, will not provide any outflow as R-H shock con¬ 
ditions are not favorable there (see Fig. 3). This certainly 
tells us that the outflow solutions would be restricted com¬ 
pared to the global shocked accretion solutions in absence of 
mass loss. However, we show that the shock induced global 
inflow-outflow solutions are not isolated solutions, but ex- 
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Table 2. Estimated Jet kinetic power from accretion rate. 


Objects 

Mbh 

(Mq) 

^acc 

i^Edd) 


£ 

ic^) 

A 

{erg) 

Xs 

(rg) 

Dmax 

^m 

(%) 

r max 
^Jet 

(ergs s“^) 

rObs # 

^Jet 

(ergs s~^) 

A0620-00 

6.60 (a) 

1.684 (b) 

0.12 (6) 

0.00444 

3.25 

26.84 

17.25 

2.48 X 10®® 



LMC X-3 

6.98 (c) 

2.487 (d) 

0.25 (e) 

0.00476 

3.15 

25.72 

17.35 

3.90 X 10®® 



XTE J1550-564 

9.10 (/) 

0.511 (g) 

0.34 (g) 

0.00530 

3.06 

22.07 

17.43 

1.05 X 10®® 

3.55 X 10®® (g) 

M33 X-7 

15.66 {h) 

0.718 (i) 

0.84t (i) 

0.00807 

2.61 

15.60 

17.79 

2.59 X 10®® 



4U 1543-47 

9.40 (fc) 

1.315 (0 

0.43 (0 

0.00565 

2.98 

21.07 

17.48 

2.80 X 10®'^ 



LMC X-1 

10.90 (m) 

0.853 (n) 

0.92t (n) 

0.00830 

2.58 

15.42 

17.93 

2.16 X 10®'^ 



Cyg X-1 

14.80 (0) 

0.061 (p) 

is 0.95! (p) 

0.00835 

2.58 

15.38 

17.98 

2.10 X 10®'^ 

3.85 X 10®'^ (y) 

Mrk 79 

5.24 X 10’^ (q) 

9.82 X 10"® (r) 

0.70 (s) 

0.00730 

2.72 

16.73 

17.68 

1.18 X lO"*"* 



M87 

3.50 X 10® (t) 

1.16 X 10-"^ (u) 

^ 0.65 (v) 

0.00705 

2.76 

16.94 

17.62 

9.25 X 10"*® 

1.00 X 10"*® (z) 

Sgr A* 

4.90 X 10® [w) 

7.89 X 10-® (x) 

0.99t (ru) 

0.00859 

2.57 

14.71 

18.09 

9.07 X 10®® 

1.00 X 10®® (zz) 

^ = Lr x Lp.dd^ where Lr is the observed Jet power (in Edd) 

(see references). 






References: (a) 

ICantrell et al\ 1 I 2 OIC 1 I) ( 6 ) IGou et al\ ll2010h (c) lO 

rosz et all ( 

2OI4I) (d) iKubota et al 

(I2OIOI) (e) ISteiner et al. 

(2014 

) 

( f) lOrosz et al 

1 (1201 la^ (oj ISteiner et alA (1201 ih (h) lOrosz et al\ 

ll2007^ fil 

Liu et all ll2^08^ (f) iLiu et all ll2ni(t) (k) lOrosz 

(200,- 

) 

(l) iMorninestar & Milled ||2014^ ( 

m) lOrosz et all ||2009^ (nj Gou et all (l2009h (0 

1 lOrosz 

et all ll2011bl) (p) IGou et all 1I2O] 

i (9) 

IPeterson et al\ 

l2004h (r) iRiffel et al\ ll2013h fsi lOallo et al\ ll201lh 

(t) 1 Walsh et all (1201,'^) (w) IKuo et all (I20l4) (ii) IWanff et al. 

0 

0 

(N 

) 


Biermann 


19991) 


ist for a wide range of inflow parameters, namely, S and A, 
respectively. Interestingly, numerical sim ulations al s o ind i- 
cate the similar findings as reported by IPas et 1 ||2014D . 
We examine the existence of such solutions with and with¬ 
out outflow and obtain the parameter space spanned by the 
£ and A as function of black hole rotation parameter au ■ As 
anticipated above, the parameter space that provides mass 
outflow is shrunk compared to the case of no outflow (see 
Fig. 7). In other words, a significant region of the param¬ 
eter space that exhibits stationary shock for Rm = 0 does 
not provide steady shock solutions in presence of outflow. 
This possibly infer that the accreting matter having £ and 
A from this region of the parameter space may demo nstrate 
the non-steady behavior (IDas et al.ll2001bl : lDasll2007h which 
is triggered simply due to the presence of mass loss. Investi¬ 
gation of such scenario requires time-dependent calculations 
which is beyond the scope of the present work and we wish 
to report this elsewhere. 

In this work, we mainly focused on thermally ultra- 
relativistic flows with adiabatic index 7 = 4/3. However, 
when the cooling effects are negligible, the purely non- 
relativistic flow behaves like gas-pressure dominated with 
7 ~ 5/3. In reality, the value of 7 would be an intermedi¬ 
ate value depending on the dissipation processes active in 
the flow. Keeping this in mind, we calculate the parameter 
space for different values of 7 and find that shock forms in 
all the cases even in presence of outflow. Interestingly, we 
observe that the effective region of the parameter space is 
reduced with the increase of 7 indicating the limited possi¬ 
bility of shock transition when the flow changes its character 
towards the non-relativistic regime. 

We have estimated the mass outflow rate using the in¬ 


flow parameters, such as, £ and A. We show that Rth in¬ 
creases with ttfe for a fixed £ and A (see Fig. 4). This is 
possibly due to the fact that as au is increased shock forms 
away from the black hole. Therefore, the effective area of the 
post-shock flow (PSC), where the inflowing matter deflects 
to generate outflow, becomes large and results enhanced out¬ 
flow rate. Further, we attempt to find the maximum mass 
loss from the disc and quantify it in terms of the 

inflow parameters. For this, we fix the value of 7 and ex¬ 
plore all possible combinations of energy and angular mo¬ 
mentum to obtain Rt^^■ In Fig. 9, we show the variation of 
with Ofe and observe that flow with 7 = 4/3 exhibits 
the highest outflow rate 77““ that lies in the range around 
~ 17% — 18%. Also, very weak correlation is seen between 
77““ and ak in all the cases. Moreover, we have shown in 
Fig. 10 that for various sets of (f, A), 77m increases with the 
decrease of Xs (equivalently Jet power increases as the size 
of the PSC is reduced) based on the realistic scenario (i.e., 
spectral states transit from LHS to HIMS) which seems to 
be a common characteristic observed in GBH sources. 

We employ our formalism in order to estimate the Jet 
kinetic power (Ljet) of several black hole sources (GBHs 
and AGNs). To begin with, we consider out bursting sources 
to calculate Lj^t based on the accretion states. For these 
sources, we compute their unabsorbed X-ray fluxes in LHS 
and HIMS using the data of RXTE observation and obtain 
their disc rate considering the accretion efficiency rj — 0.3 
which seems to be relevant for rotating black holes. Then, we 
calculate the outflow rate using our formalism and employ 
it in Eq. (15) along with the disc rate to obtain Lj^t- While 
estimating the Jet kinetic power, we consider maximum out¬ 
flow rate 77)0“ for HIMS (see Fig. 9) and 77m ~ 10% for 
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LHS, as discussed in Section 5. In Table 1, we summarize 
the physical parameters of five black hole sources along with 
the computed Jet kinetic power. In the process of estimat¬ 
ing UjNj f"'® quantities play major role. First one is the 
unabsorbed X-ray flux computed from RXTE archival data 
and the other is the maximum outflow rate obtained from 
our model calculation. We find that the estimated UjN for 
various sources are in clo s e agreement with the observed 
values l|Fender et al.l[2003 : iMiller et al.ll20f3 l. We further 
extend our study for other black hole sources where the fast 
accretion state transitions (~ day scale) are not seen (ex¬ 
cept XTE J1550-564). For these sources, the disc rate is ob¬ 
tained from the existing literature (see Table 2) and 
is computed from our theoretical calculation which finally 
provides the We present these results in Table 2. We 

notice that estimated UjN for Cyg X-1, XTE J1550-564, 
M87 and Sgr A* are also in close agreement with the ob¬ 
served values (iF ender et al.l l200^ : Ide Gasoerin et ^l2012l : 
iFalcke fc BiermanrJ 1999l l. Following this findings, we argue 


the representative values which are expected to be consistent 
with their actual estimates. 

The present work has limitations as it is developed 
based on some approximations. We use pseudo-Kerr grav¬ 
itational potential to mimic the general relativistic effect 
around rotating black hole that allows us to examine the 
properties of non-linear shock solutions in presence of mass 
loss in a simpler way. We ignore the effect of viscosity and 
radiative processes. We consider constant adiabatic index in¬ 
stead of calculating it self-consistently based on its thermal 
properties. Needless to mention that we have not address 
the issue of transient relativistic ejections and its collima- 
tion mechanism as we consider the Jet geometry only up 
to its sonic point. Although the implementation of all such 
issues are beyond the scope of this paper, however, we be¬ 
lieve that the above approximations will not alter our basic 
conclusions qualitatively. 
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